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The four- dimensional A^A*" — nNN equations are adapted to the case of scalar particles with 
a 0^(T interaction Lagrangian, and solved for energies below the (T-production threshold. This is 
achieved in the approximation where <j!>(t scattering is dominated by the s-channel (/>-pole term. 
The importance of the removal of double-counting is investigated, and a detailed comparison of the 
results of a covariant coupled-channels formulation and the Bethe-Salpeter equation in the ladder 
and ladder plus crossed-box approximations is presented. A brief discussion of the extension of the 
method to energies above the cr-production threshold is given. 



I. INTRODUCTION 



Quantum Chromodynamics (QCD) is now the accepted theory of the strong interaction. Its highly non-perturbative 
character at hadronic mass scales means that effective field theories of the strong interaction are almost always 
■ used when the properties of few-hadron systems are calculated. In some such field theories (e.g. those based on 
■r^-j- ' Chiral Perturbation Theory) the Lagrangian is organized so that a perturbative calculation of hadronic reactions 
is appropriate, at least for "small" momenta [Q-^. However, to examine the scale at which quark-gluon degrees 
I ' of freedom become essential we need to use these field theories to calculate hadronic reactions at medium to large 
I momenta. The need for explicit quark-gluon degrees of freedom may then be deduced from failures of effective 
' hadronic field theories to explain the experimental data. It is in this regime that Chiral Perturbation Theory fails, 
and unitarity plays an important role. 

The amplitudes for any hadronic reaction are, of course, exactly given by the infinite hierarchy of Schwinger-Dyson 
' equations, which can be derived by functional techniques An alternative method of deriving a set of coupled 

\^ . equations for the amplitudes of a field theory was developed by Taylor j^J^ . It involves taking the series of Feynman 
On ■ diagrams for a particular amplitude and classifying each graph according to its topology. The result is an infinite 
, set of coupled equations for the amplitudes of the field theory. Some approximate truncation of this set of equations 
• (e.g. one which preserves s-channel unitarity) must be found. Since these equations are derived directly from the 
field theory, all of the integrals present in them will be four-dimensional. Note that despite the use of the Feynman 
diagrammatic series, Taylor's method is a non-perturbative technique, since it sums certain classes of graphs to all 
orders in the coupling constant. The great advantage of this method is that the resulting equations depend only on 
the topology of the allowed vertices of the field theory. Hence this is a particularly useful approach in hadronic field 
theories, where the best effective Lagrangian has not yet been derived from QCD. 
, When applied to the two-to-two amplitude the Taylor method yields 

tW =r(2)-f r(2)did2T^i), (i.i) 

where the bracketed superscripts indicate the irreducibility of each amplitude. This is, of course, a Bethe-Salpeter 
equation (BSE) for the two-to-two amplitude, with T^^^ as the "driving term" We note that if T^^^ is the sum 

of all two-particle irreducible two-to-two Feynman diagrams, then this equation is exact. If the Taylor method is 
applied to the three-to-three amplitude, four-dimensional Faddeev-like equations may be derived |11|. 



In hadronic physics one problem which one might try to apply Eq. (1.1) to is nucleon-nucleon scattering. Before any 
calculations can be performed the nature of the kernel T^^^ must be specified. One approach, pursued by Tjon and 
Fleischer, is to assume T^^^ consists of a sum of one-boson-exchange Feynman diagrams. This approximation gives a 
reasonable fit to the nucleon-nucleon phase shifts ||l^,|l^. However, for energies above the pion-production threshold, 
one needs to include the thresholds for all of the physically allowed channels. In particular, if a theory which gives 
consistent predictions for NN scattering and pion production is to be developed, the full NNn intermediate-state 
structure of the amplitude T^^^ must be exposed. The Taylor method is well suited to this task, and attempts to 
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use it in order to expose the NNtt intermediate-state structure of field-theoretic amplitudes, such as T^'^\ were first 
pursued by Avishai and Mizutani jlj,|l5|. They used Taylor's method to classify all of the covariant perturbation 
theory diagrams which contribute to the NN NN, NN ^ NNn and NNtt NNn amplitudes and contain one 
explicit pion. The result was a set of Faddeev-like coupled four-dimensional integral equations for the processes 

N+N\ ( N+N , . 

N + N + 7:j^\N + N + n ' ^ 

These equations were formally equivalent to the three-dimensional NN — ttNN equations previously obtained by a 
number of groups using different techniques However, a crucial difference was that Avishai and Mizutani's 

equations contained four-dimensional, rather than three-dimensional, integrals. Similar four-dimensional equations, 
but for the BB — ttBB system, where the B could be either a nucleon or a delta, were later obtained by Afnan 
and Blankleider |23| . By explicit construction, both Avishai and Mizutani's and Afnan and Blankleider's equations 
were covariant. Furthermore, unlike the three-dimensional NN — ttNN equations, they did not treat different time- 
orders of the same physical process differently. Consequently, these four-dimensional equations avoided the theoretical 
problems of the three-dimensional NN — ttNN equations which were pointed out by Sauer, Sawicki & Furui and 
Jennings & Rinat [^,|6). (For a more detailed discussion of this point see |27|j2^ .) 

Unfortunately, these four-dimensional NN — ttNN equations were not correct, since they double-counted certain 
Feynman diagrams. This double-counting occurred due to a problem with the original Taylor method. Indeed, 
such double-counting will arise whenever one tries to derive coupled field-theoretic equations whose input is full sub- 
amplitudes for sub-system interactions. While a decomposition in terms of sub-system interactions is possible for 
any Feynman diagram, it is not necessarily unique. Therefore, when the perturbation series is reorganized in order 
to derive such coupled equations, some Feynman diagrams get incorrectly included more than once in the resummed 
series. (See |29j for a brief discussion and ||] for details.) To overcome this problem, we need a procedure for 
determining the conditions under which double-counting occurs, and a way to remove the offending diagrams from 
the equations. A "modified Taylor method" , which gives a prescription for the subtractions necessary to remove 
this double-counting, was recently developed by us Q. We then used this modified Taylor method to derive revised 
four-dimensional NN — ttNN equations which do not contain double-counting [^,p9| . These equations are essentially 
equivalent to those obtained by Kvinikhidze and Blankleider using different means | |3C| , |3l| ] . The resulting set of coupled 
four-dimensional two-fragment scattering equations obey NN and NNtt unitarity, are free of double-counting, and 
are covariant, not only on-shell, but also off-shell in the manner dictated by the Feynman diagrammatic expansion. 

Thus a set of equations which gives a complete and correct summation of the NN and NNtt sectors of any hadronic 
field theory now exists. However, it remains to be seen whether this set of equations can be solved in order to yield 
predictions for experimental quantities. For a detailed program describing how such a solution might be achieved, see 
Section 7 of Ref. [ p7| . This paper reports on a preliminary study of the numerical solution of these equations. The 
equations are solved, not in a field theory of nucleons and pions, but in a scalar field theory. Consequently, details 
of spin and isospin are removed from the problem. In addition, there is no need to introduce subtractions or form 
factors to obtain convergent integrals. The simplified calculations described here imply that the numerical solution 
of the full four-dimensional NN — ttNN equations is feasible. The use of a scalar field theory also makes for a clear 
discussion of issues such as the usefulness of a coupled-channels formulation and the importance of double-counting 
removal in such an approach. 

A field theory describing two types of spin-0 particles interacting via a Yukawa coupling is chosen. A contact 
interaction is also included, in order to ensure the stability of the classical vacuum. The Lagrangian density for this 
"<^V field theory" is 

£ = liid.c^r - m'c^' + {d,af - ^?a^] - f^V - ^^^a^ (1.3) 

Because the modified Taylor method is Lagrangian independent, the NN — ttNN equations for distinguishable 
particles derived in Ref. |^ apply to the (fP'a field theory. The distinguishable-particle equations must be symmetrized, 
as the (/)'s are bosons. This leads to a set of nine two-fragment scattering equations for the processes 

(j) + (t)\ [ (1) + ^ 

{cj^a) + ^{ (0a) + , (1.4) 
(00)d + cr J [ {(j)(f))d + (7 

where {4>4>)d is a bound state of the (j)(j> system that plays the role of the deuteron. Below the threshold for a 
production these equations provide a way of calculating the </> — </) amplitude which corresponds to including infinitely 
many diagrams in the kernel of a single Bethe-Salpeter equation for that amplitude. Above cr-production threshold 
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these coupled equations allow us to calculate amplitudes for the reactions of Eq. (1.4) in a framework which — provided 
the input amplitudes have two-body unitarity — maintains (jjcj) and 4>(t)a unitarity. Note that as in calculations using a 
single BSE, the kernel to these equations must be approximated in order to make them computationally tractable. 
The inputs to these — acfxj) equations are the one-particle irreducible a — <j) amplitude t^^J and the (j)(j)a vertex 

function /'•^■'. Because the (fxf) — > (fxf) amplitude appears in the kernel of these equations they are a set of non-linear, or 
bootstrap, equations. Strictly speaking the complete non-linear set of four-dimensional equations should be solved for 
the amplitudes of interest. To snap the bootstrap, in this work we mimic the approach of previous three-dimensional 
■kN — ttttN calculations and use an input (fxj) amplitude which is constructed as a separable approximation to the 
"true" (p ^ 4> amplitude. The energy range used in the construction of this expansion must be /i below the energy 
at which the coupled equations are to be solved. Arguments from Ref. ||2^ show that we expect such an approach 
to yield a reasonable solution to the set of non-linear coupled equations. In this work we take the predictions of the 
ladder Bethe-Salpeter equation as the "true" — </> amplitude. The consistency of this procedure may be judged by 
comparing the result of the coupled-channels calculation to the separable input amplitude. 

(2) 

In order to simplify matters as far as possible we assume that the amplitude r^^ is zero, i.e. that tpa scattering is 
dominated by the pole diagram. The implications of imposing this restriction are discussed in Sec. II. This assumption 
implies that the vertex function for a absorption on a single (p is given by 

f^'\p\p,k) = f^^\p',p,k) = -ig{2n)H('\p' -p-k), (1.5) 

where the Feynman rules of the field theory have been used in order to evaluate /(^^. 

If the input amplitudes are chosen in this way then the (jxj) — a(j)(j) equations become a set of four coupled equations 
describing the processes in ( |l.4D which do not involve the (0a) "isobar" . In this paper we solve these equations for 
(f>4> scattering up to the cr- "deuteron" threshold and compare the results with those from the BSE. In so doing we lay 
the foundation for studies of the covariant, four-dimensional, NN — ttNN equations derived in Ref. |^7|] . 

In Sec. II we consider the full (j)4> — <j4>(j) equations and discuss the consequences of the approximations made to 
get the simplest form of the equations which includes the coupling to the a{4'(j))d channel. In the process we see how 
the double-counting problem arises at each level of approximation, and how the content of the kernel compares with 
that of the BSE. We also examine the consequences of taking the ladder approximation for the 4> — (j) amplitude in 
the kernel of the coupled integral equations. 

Since the input to the equations derived in Sec. II is the (j) — (j) amplitude, we proceed in Sec. Ill to discuss the 
numerical solution of the Bethe-Salpeter equation in the ladder approximation. We demonstrate that, in the case 
when the interaction supports a bound system with a binding energy comparable to that of the deuteron, a valid 
separable approximation to the — </> amplitude can be written. 

In Sec. IV we examine the solution of the coupled integral equations for the reactions in Eq. (1.4). At this stage 
only the solution of these equations below the threshold for a production is explored. We then compare our results 
with those from the BSE when the kernel is taken to be either the ladder, or the ladder plus crossed-box a exchange. 
The effect of vertex dressing in the BSE is also investigated. In particular, it is found that vertex dressing must be 
included in the Bethe-Salpeter kernel in order to maintain consistency with the coupled-channels results. Finally, in 
Sec. V we give a summary with some concluding remarks about future work. 



II. THE (t)(j) - (T<i)(l> EQUATIONS 

There are two central questions in the analysis of the ac/xj) equations. Firstly, what additional features do these 
equations have as compared to BSE formulations? In particular, how does the physics content of these equations 
change as different approximations are made for the kernel? Secondly, what numerical problems are encountered in 
solving these equations? In this section we address the first of these questions, by illustrating the origin of the double- 
counting problems and showing how subtractions remove them. The NN — ttNN equations reported in Ref. |^,^ , 
which are basically the same as the equations of Kvinikhidze and Blankleider [popi|] , are also valid for the Lagrangian 



in Eq. (1.2) since the derivation of the equations does not depend on the choice of Lagrangian. As a result, the 



acfxf) equations for distinguishable (/)'s are of the form 

= V {1 + did2 T^4, ) 
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a 

- h f^^^Hj) ( 1 + did2 T^4, ) , (2.2) 

where we have taken particles one and two to be the (j) particles, while particle three is the a. Here, 6ij = 1 — 6ij and 
j — 2 for j = 1, while j — I for j ~ 2. The roman indices run from one to two, while the greek indices run from 
one to three. These equations differ from the corresponding equations of Avishai and Mizutani MM and Afna n and 



Blankleider [gsj in that the (f>—4> potential V and the one-particle irreducible 2^2 amplitude i'^^^^iafin Eq. (2.1) have 
subtractions to remove double counting in the (f> — (j) amplitude T^^, while the last term in Eq. ( |2.2| ) is a subtraction 
to avoid double counting in the cjjcj) acfxp amplitude T\^. /'■^•'(i) is the one-particle irreducible (p 'r— (f>a amplitude 
with particle i as spectator, while t'^'(a) is the one-particle irreducible 2^2 amplitude with par ticle a as spectator. 



The single-particle Feynman propagator for particle a is do,. The a — (f> potential in Eq. (2^) corresponds to the 
M-channel 0-pole diagram, i.e. 

^;-(z) = /«(z)d, /«t(^) . (2.3) 

The two-body amplitude i'^^^ is chosen such that the subtraction is present only in the a — (f> amplitudes (see Ref pTf), 
with the subtraction basically removing any initial u-channel pole, i.e. 

^2=^i; i':h{t'-ld,d. + l)v^,; (2.4) 

«?0=iS-<0- (2.5) 

In Eq. ( 2A ) the 4> — (j) potential V is given by 

y = V5sE -Di-D2-D,-X-B , (2.6) 

where Vq^^ is the one- cr- exchange potential, i.e., 

VSsE^f'-^^*{^)daf^^^*\2), (2.7) 

with the modified one-particle irreducible cri^^ form factor, f^^^*, having dressing with the w-channel pole diagram 
subtracted, i.e. 

+ ■ (2-8) 



The subtractions D2, Da-, X, and B are illustrated in Fig. and are given explicitly in Ref |27 . Di is merely D2 

with the two nucleons interchanged. All sub-amplitudes in Fig. |l| are one-particle irreducible in the s-channel. 

Since we are ultimately interested in the NN — ttNN equations, the first approximation we consider is that the 

(2) 

two- particle irreducible cr — amplitude is just the u-channel pole diagram, i.e. r^! = v^i. In this approximation 



the subtractions Di and D2 in (^) become unnecessary. All of the other double-counting problems are still present, 
and must be removed by the explicit subtraction terms. Since v^^ = 0, the form factors in the cr-exchange potential 
are two-particle irreducible, while the u-channel subtracted a — (j) amplitude is zero, i.e. 

We could now approximate both the a — (j) and the (p — (j) amplitudes by separable expansions, which, after partial- 
wave expansion, would reduce the coupled-channels problem to a set of coupled integral equations in two dimensions. 
However, since this is the first numerical work on this problem, we have chosen to concentrate on the problem of the 

separable expansion for the 4> — (f) amplitude, and the procedure for linearization of the coupled non-linear integral 

(2) 

equations. To this end we have chosen t 1 = 0. Such a choice implies that all one-particle irreducible form factors 



in Eqs. ( |2.l|) and (2^) be replaced by the corresponding two-particle irreducible form factor, and the one-particle 

irreducible a — (j) amplitude t^^^ be set to zero, i.e., f^^'^ /'^' and t'^-* —t 0. This reduces our coupled integral 
equations to two coupled equations of the form 
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T4>4> = V ( 1 + did2 T^4, ) 

+ E /'"^ (j') '^J^- *S ^1^2 Trf^ (2.10) 

Td^ = ^ /^'^^(j) ( 1 + did2 ) . (2.11) 

In writing these equations we have made use of the notation Td^ = Ts^. In these equations the correction required to 
remove the double-counting problem is restricted to the (j) — (f) potential V . This is because all double counting due 
to u-channel poles in the a(j) t-matrix has been, by definition, removed from the problem. 

At this stage we should recall that the one-particle irreducible (/) — 4> amplitude, t^^^ , is equal to the full amplitude, 
T^^, which is part of the solution to the coupled equations, and so we have a non-linear set of integral equations. 

In the time-ordered, or non-relativistic, version of these equations, the energy variable in t^^^ is lower than the 
energy in the amplitude T^^ by at least the mass of the a. This allows for a bootstrap procedure, in which the 
lower energy amplitude is used to obtain the amplitude at higher energies — as illustrated in work on the nN — mrN 
equations above the threshold for pion production . A similar situation prevails in the covariant framework — but 
only approximately (see ||2^ , ^^ for details). In the present analysis we will follow this procedure for removing the 
non-linearity. To simplify the resulting bootstrap problem we will assume that t^^^ is a solution of the BSE in the 
ladder approximation. This final approximation removes the bootstrap problem, and our coupled-channel equations 
still satisfy two and three-body unitarity. Furthermore, we can test the validity of this solution to the bootstrap 
problem by comparing the phase shifts corresponding to i^^j with the phase shifts resulting from the solution of the 
coupled-channels problem. 

Before we turn to further discussion of the input to the coupled-channel equations, we should point out that the 
above two coupled equations can be recast into a single equation for the 4> — (f) amplitude 

= K// ( 1 + did2 T^^ ) , (2.12) 

with an effective interaction of the form 

Veff = ^ + E /^'^('^) did. t^ll f'-'^Hj) ■ (2.13) 

Rewriting the coupled-channels problem in terms of a single channel allows us to see how the double counting is 
removed by subtractions in V. It also helps in discussions of the difference between the coupled-channel equations 
and the ladder, or ladder plus crossed-box, BSE. 

In the event that t^i chosen to be the solution of the ladder BSE, then the second term on the right hand side 



(r.h.s.) of Eq. ( ^.13 ) has two terms that correspond to crossed-box a-exchange. This double counting is removed by 
taking 

V = VSsE-B , (2.14) 



where B is the crossed-box digram (see Fig. |i|). In addition, when i = j, the second term on the r.h.s. of Eq. ( 2.15 ) 



has diagrams that contribute to vertex dressing of the cr-exchange potential Vq^^ (diagrams (a) and (6) of Fig~|) 
However, a careful examination of the series resulting from this term reveals that there are no contributions from 
diagrams in which both vertices are being dressed at the same time. This is a result of our elimination of the non- 
linearity in the coupled integral equations. Finally, expanding the i ^ j piece of the second term on the r.h.s. of 
Eq. ( p.l3| ) shows that some of the cr ossed -ladder diagrams required to get the correct one-body limit for the BSE |3~ 



are included in the solution of Eq. (2.12). Indeed, the solution of the non-linear equation for T*^^^ will include more 
(but not all) of the crossed-ladder diagrams necessary for the corre ct one-body limit. 
In the event that the amplitude from the solution of Eq. ( 2.12| ) is not consi stent with the solution of the ladder 



BSE, we could take the solution of Eq. ( 2. IS ) and su bstit ute that result in Eq. ( 2.13 ). This in turn will give us a new 



effective potential that can be use d in the BSE, Eq. (2.12). However, in that case more diagrams are generated by the 
second term on the r.h.s. of Eq. (|2.13| ), and therefore further subtractions are needed. . Consequently, the potential 
V is given by 

y = y5sE -B-D,-X . (2.15) 
This demonstrates very clearly that when the input to the equations is approximated we need to be careful in 



calculating the correction for the double-counting problem. We return to this point in Sec. IV D when we consider 
including vertex dressing in the BSE. 
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III. THE BETHE-SALPETER EQUATION 



In most numerical analyses of the NN — t:NN equations the input two-body amplitudes were parameterized to 
fit the experimental data, and no regard was given to the consistency between the input N — N amplitude and the 
N — N amplitude resulting from the solution of the coupled-channels problem. In the present analysis there is no 
experimental data for the 4> — 4> system, and we would like to specifically maintain some consistency between the input 
4> — (ji amplitude and that resulting from the solution of the coupled-channels problem. Therefore, as we have detailed 
in the previous section, we will consider the coupled <f)(l) — a(j)(j) equations with the only input amplitude being the (j^—cj) 
amplitude in the ladder approximation. In this way the problem of the non-linearity of the coupled-channel equations 
is avoided. After partial-wave expansion, the resultant coupled-channels problem is a set of coupled integral equations 
in four dimensions. To reduce the dimensionality of our equations we follow the procedure previously implemented 
for the TriV — ttttN equations |^^, and construct a separable approximation to the solution of the ladder BSE. In this 
way we reduce the dimensionality of the coupled integral equations from four to two, at the possible cost of increasing 
the number of coupled integral equations. 

Therefore in this section we explain various aspects of the solution of the Bethe-Salpeter equation and the construc- 
tion of our two-body input. These include techniques for: the solution of the ladder BSE; calculation of bound-state 
masses; construction of a separable approximation to the ladder amplitude; and inclusion of effects beyond the ladder 
approximation. While all of the techniques described here are explained elsewhere ^J-Q, we give an outline of the 
methods for two reasons. Firstly, our solution procedures for the coupled-channel equations are extensions of those 
described here. Secondly, explaining the methods establishes some of the notation used in Section IV. Further details 
on the work of this section may be found in Ref. p8[ . 



A. The ladder approximation 



The BSE, Eq. (1.1), is an equation for the invariant two-particle scattering amplitude in a field theory. Of course, 
in practice it is impossible to construct the sum of all graphs contributing to T*^^^ and solve the resulting integral 
equation for T^^\ Therefore, most work in the past has focused on the so-called "ladder approximation", in which 
T^^' is taken to consist solely of one-particle exchange. 

Consider the (fP'a field theory. Define: 

T^^\p\,p'^■,Pl.P2)^-^{2^rf5(^\p',+p',_-p,~P2)T{<^,q■P) . (3.1) 

If the total momentum is P = p'^ + P2 = Pi + P2, then in the kinematics specified by Fig. || the BSE in ladder 
approximation is 



r(g', q- P) = Viq', q) + J d^q" V{q', q") d^{ql) d^{q'L)T{q'\ q; P) . (3.2) 

Here = ^ ± q" , while the potential due to a-exchange and the propagator are defined by 
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^ii''i) = -ri — ^^2 — 2' Mp) = ^-^ — 2- (3-3) 

[q' — qY — fi'' p^ — 

Throughout this paper the masses m and p are regarded as having an infinitesimal negative imaginary part. This 
defines the way poles should be negotiated, unless otherwise stated in the text. 

The BSE may also be used to extract information about bound-state properties. It is known that if a bound state 
of the two (f) particles with mass M exists then the invariant scattering amplitude may be decomposed into a pole 



and non-pole piece Inserting this decomposition into Eq. (O) and taking the residue at the pole yields 



r(g';M') d'q" V{q\q") [d^{ql)d^{q'L)\p^.^,,M';M^) , (3.4) 

where F is the bound-state vertex function. Thus the bound-state solutions of the ladder BSE may be found by 
searching for those values of P^ for which the kernel of the BSE in the ladder approximation has an eigenvalue of one. 

The BSE in the ladder approximation is a singular integral equation and as it stands is not amenable to numerical 
solution. One possible way of dealing with the singular nature of the equa tion is to perform a Wick rotation pM , which 



involves an analytic continuation of the two variables q'^ and q^ in Eq. (3.2) to the imaginary axis. In the center-of- 



mass frame, where P = {\fs, 0), this analytic continuation requires that we take into consideration the singularities of 
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both the kernel and the Bethe-Salpeter ampHtude. The singularities encountered depend on the total energy available, 
y/s. As explained by Taylor and Paganamenta this straightforward analytic continuation to the imaginary axis 
is valid for < < 2m + 2^. Above -/s — 2m + 2/i the self-consistently-generated (j)a production threshold cuts 
in the amplitude T, which begin at Qq = ±(-^ — m — ji), enter the first and third quadrants, thus pinching the Qq 



integration contour, and preventing a simple analytic continuation of Eq. (3^). Note that for 2m < ^/s < 2m + 2/x 
poles of the (f> propagators may enter the first and third quadrants of the q'^ plane. If this occurs then the residue at 
these poles must be included in the Wick- rotated equations, in order to ensure the correct analytic continuation. 



To reduce the dimensionality of Eq. (3.2) from four to two, and so be able to turn the resultant equation into an 
algebraic equation, we need to partial-wave expand both the amplitude T and the one-cr-exchange potential. The 
partial-wave expansion 

A(g;,g^;s) = (-Stt^)^ A,{q'o,q',qo,q;s) Pi{cos0) , (3.5) 

where A is either the one-cr-exchange potential V or the amplitude T, and cos 6' = q-q', is used. In the present analysis 
we take fi = c = 1, and the unit of energy is chosen such that m = 1. From this point all four- vectors have a subscript 



fi, so as to distinguish them from the magnitude of three- vectors. The partial-wave expansion converts Eq. (3.2) to 

oo oo 

Te{q'o,<l',<lo,q;s) ^Vi{q'Q,q',qo,q)-i J dq'g [ dq" Vi{q'„, q' , q'^ , q") G{qo , q"; s)Ti{qa, q" , qo, q; s) , (3.6) 

— oo 

where: 







Gilo, q": s) = y= , (3.7) 

(^ + q'^)^ - E^q")' (f - q'^r - E^{q"Y 

' ^ ( q'^ + q^+li''-{q'a-qo? \ ,„ 

TT V 2q'q J 

E^iq")^y/q"^+m^~ie; A = ; (3.9) 

and Qe is a Legendre function of the second kind. 

The on-shell partial-wave amplitude T|'" can now be written in terms of the phase shift di as 

Tr{s)=Te{0,q,0,q;s) = ^ e^^^(^) smS,{s) , (3.10) 



where the on-shell relative momentum is given hy q — — tiv^- 

Next we perform a Kowalski-Noyes (KN) p^-p^ reduction, in order to remove the two-body unitarity cut, i.e. we 
write 

T,{,^^,q" ,0,q;s) = fi{q'^,q"-s)Tr{s) , (3.11) 

with the KN half-off-shell function being one at the on-shell point, i.e. fi{0,q;s) = 1. This is followed by a Wick 
rotation. Since the two (f> particles have the same mass, we change the q^ integration to one from zero to infinity, and 
so obtain the coupled equations derived by Levine et al. for fi{iqQ, q'; s) and ge{q'; s) p^ , ^ . The auxiliary equation 
for ge{q'', s) must be written because of the presence of 0-propagator poles in the first and third quadrants of the q'g 
plane. gi{q; s) is defined by 

gdq; s) = Mw{q), g; s) , w{q) = — ~ E^{q) . (3.12) 

The equations for fi and gi can be solved by discretizing the integrals and using direct matrix inversion. We note 
that there is a term Vi{w{q'), q', —w{q"), q") in the kernel of the integral equations. For 2m -I- /i < y/s < 2m -I- 2/i this 
term may have a logarithmic branch-point in the region of integration. We found that the equations could still be 
solved accurately in this energy regime, provided enough quadratures were used for the discretization. 
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To calculate phase shifts, the equation 



' ^ ' 1 + * /"oo /o°° ^9" ^^(0, 9" g^,', 9") G(g^,', 9"; s) Mq'^, q"; s) ' ^ ' ' 

is used to obtain T^°^^(s). In order to deal with the singularities in the integrand the integral is written as 

oo oo 

i [ dq'^ [ dq" VfXO, q, q'^, q") G{q'^, q"; s) M, q" ■ s) = 2a{s) + h{s) , (3.14) 



where the factors a{s) and h{s) are given by 



oo oo 

a{s) - - I dq'^ j dq" G{iql q"; s) (T/,(0, q, iq^ q") h{iql s) - Vr{s)) 


+ TT I'dq" ] {V,{Q,qMci").q")9i{<l"\s)~Vr{s)) , (3.15) 

Jo 2yf^w{q")E^{q") 

OO OO 

b{s) = Vris) j dq'^J dq" G{q'^,q"-s) . (3.16) 

-oo 

To test the accuracy of our numerical procedure, we compared our results with those reported by Levine et al. 
and found agreement. Varying the quadrature distribution then shows that with 14 and 22 Gauss-Legendre quadra- 
tures for the (7o and q" integration in the Wick rotated integral, and 60 in the auxiliary residue q" integral, the real 
part of the phase shift is accurate to three significant figures, while the imaginary part is accurate to two significant 
figures. The large number of quadratures is necessary in the auxiliary integral because of the presence of logarithmic 
branch-points in the kernel if -^/s > 2m -I- /i. 

Since we hope to have in the 4)(t> — a(l)(j) system a simple model for the NN — ttNN system, we have chosen to 
construct the input (f) — (f> amplitude to have a bound state with a binding comparable to that of the deuteron. To 
achieve this we consider the homogeneous BSE in ladder approximation after partial-wave expansion, 

oo oo 

Ti{q'o,q';s) ^-i I dq'^ [ dq" V,{q'o, q' , q^ q") G^q^q" ; s)T,{qlq" ; s) , (3.17) 



where Ti^q'^, g'; s) is the vertex function for relative four-momentum g^^, with its angular dependence given by Pi{q' -z). 

This may be rewritten as an eigenvalue equation with a symmetric kernel, Ki. Note that because is linear in the 
coupling A it is not actually necessary to calculate the eigenvalue at a number of energies. Instead we merely write 
Ki = X K[. The largest eigenvalue of K'^{s) is then the inverse of the coupling A required to produce a bound state of 
energy s with angular momentum £. Using this technique we obtained, for the case /i = 0.15, the plot of bound-state 
energy against A shown in Fig. IJ. If a result accurate to three decimal places is desired the q'g and q" integrations 
require 12 and 22 quadratures respectively. 

The coupling A is chosen so that the S'-wave 4> — (j) bound state has the same binding energy as the deuteron. In 
our units = 3.9905. The value A = 0.131 yields a bound state of this mass. The value of the "deuteron" binding 
energy is only accurate to 5%; however, for our purposes this is more than adequate. 



B. Separable approximation 

Having determined the Bethe-Salpeter amplitude in the ladder approximation, we now need to get a separable 
representation of this amplitude, so that we can solve the (jycf) — a(j)(f> equations numerically. Ideally, a separable 
expansion of the one-cr-exchange amplitude in terms of the eigenstates of the kernel, or the half-off-mass shell solution 
to the BSE in the ladder approximation |^], should be used. However, as a first attempt at a solution of the 
— a(f>(j) equations, we consider the phase shifts from the ladder BSE as data, and fit this data with a rank-one 
separable potential. 

In a given partial wave, we can write a covariant version of a rank-one separable potential, Vi, as [56| 
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Vtiq'o, 90, q) = -^qxiiq'o^ q')it xeiqo, q)q ■ (3.18) 



The factors of q and q' are introduced into the potential to maintain consistency with our BSE, Eq. (|3.6|) . Since we 
will be restricting our analysis to S'-wave scattering, we have chosen the form factor for ^ = to be a function of 
the four-momentum square, 

Xo(go,9) = 2 _ \ — M ■ (^-i^) 
90 9 P 



Te{q'o,q\qo,q;s) = --^^q' Xtiq'o^q) Te{s) Xi{qo,q) q , (3.20) 



The solution of the the BSE for this potential is also separable, being of the form 

1 

47^3 

where: 

oo oo 

r,-\s) = ^ - / ^90 / dq"[q"xi{q'^, q")? G{q'^, q"; s) . (3.21) 

-oo 

Note that, assuming the existence of only one bound state in the channel with angular momentum £, ti may be 
rewritten 

Ms) = ^ , (3.22) 
where is the bound-state mass, t o be d etermined by demanding that ti{s) has a pole at s = m^, and Se{s) is the 



residue of ti(s) at that point. In Eq. ( 3.21 ) the integration is performed analytically using the residue theorem 
For the scattering case, to remove the pinch at qQ = 0, q" = q, we rotate the contour of integration into the third 
quadrant of the (/"-plane. An angle of rotation of ^ and 32 quadratures for the q" integration gives sufficient accuracy 
to fit the phase shifts resulting from the solution of the ladder BSE. 

The parameter is adjusted to guarantee that the binding energy is identical to that predicted by the BSE 
equation in the ladder approximation with /x = 0.15 and A = 0.13, while the parameter f3 is adjusted to get a best 
fit to the phase shifts below the production threshold. The final values for the two parameters in our choice of units 
are: = —15.33 and f3 = 0.47386. Despite our use of only two parameters, the fit to the 'data' is very good, as 
demonstrated in Fig. |[ 



C. Beyond the ladder approximation 

In this section we consider two improvements to the standard ladder approximation, which will be used when we 
are comparing single-BSE calculations with results from the coupled-channels calculation. Firstly, we examine the 
inclusion of (j) propagator dressing, which is necessary if cjycfxj unitarity is not to be violated. Secondly, we explain how 
to include crossed-box cr-exchange in the kernel of the BSE. 

The first improvement is to include the lowest-order dressing for the (j) propagators. This effect was included 
in the derivation of the (j)4' ~ '^4'4' equations by implicitly assuming all propagators were dressed. In addition, the 
dressing of the propagators gives a contribution to three-body unitarity that complements the contribution of one-cr 
exchange. Indeed, Levine et al. have shown that such propagator dressing must be included in order to guarantee 
an inelasticity less than one above the threshold for a production. They found that for fi = m = 1, and sufficiently 
large coupling, the dressing of the 4> propagator plays an important role. 

The inclusion of the minimal dressing in the propagator takes the form (see Fig. |^) 

<P^) = 2 2^ ^ ■ ' (3-23) 



where to lowest order in the coupling 

d^q 



^ (p — qY — /i^ q^ 



(3.24) 
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Note that we have not included the full one-loop self-energy term here, since in order to satisfy (j)<f>a unitarity it is only 
necessary to use bare propagators in the expression for E. We now write a dispersion relation for this propagator. 
Since we wish to examine the propagator after mass and coupling-constant renormalization this is a twice-subtracted 
dispersion relation. Algebraic manipulations a la Saenger |^ and Janus ||5^] yield: 

'^^P^ ~ p2 -m? + i€ l + [p^ - m?)A{p'^) ' (3.25) 

with: 

oo 

A{p)^\ / ds ^ — — 2 n • (3-26) 



(m+/j)2 



The integral (3.26) may be evaluated numerically. We find that 18 quadratures are more than sufhcient to ensure 
that it is accurately calculated. 

Previously, Levine et al. have shown that for m = fi — 1 and a large coupling constant, the effect of the 
dressing on the phase shifts is substantial. Although we agree with their numerical results, we find that for m = 1 
and /i = 0.15, which correspond to the mass of the nucleon and pion, and a coupling strength that gives a binding 
energy for the — </> state that is comparable to the deuteron binding, the effect of the dressing on the phase shifts is 
small. This is particularly the case for energies below the tr-production threshold, as demonstrated in Fig. ^, where we 
present the real part of the phase shifts. In fact, with dressed propagators a coupling of A = 0.13 yields a bound-state 
of mass — a very small change from the result for undressed propagators. 

Despite this dressing having very little effect on the bound-state position and the real part of the phase shifts, it does 
have a substantial impact on the imaginary part of the phase shifts, as seen in Fig. ^j. However, the imaginary part of 
the phase shifts is still small. This result is consistent with observations in N — N scattering above pion-production 
threshold, where the contribution to pion production is dominated by the production of a A which then decays by 
pion emission, and the equivalent mechanism via the nucleon is comparatively small. 

A more interesting addition to the ladder approximation is the inclusion of the crossed-box cr-exchange diagram 
in the potential. In particular, the coupled-channel equations with the ladder approximation for the input (p — 4> 
amplitude include such crossed diagrams, as pointed out in Sec. |^ Furthermore, the inclusion of the crossed-box 
graph in the kernel partially corrects the failure of the ladder BSE in the one-body limit. If the BSE is written for 
two particles with unequal masses, then when the mass of one particle is taken to infinity we expect to recover an 
equation for the other, lighter, particle which represents it moving in the static potential generated by the heavy 
particle. If the ladder BSE is used this does not occur. This one-body limit is restored on including all ladders and 
crossed- ladders in the perturbation series for the amplitude Although the inclusion of the crossed diagram in 

the potential is not equivalent to including all crossed-ladder diagrams, the fact that our coupling constant is small 
implies that including this one extra diagram in the kernel will give us an indication of the contribution of higher-order 
crossed-ladder graphs. 

Adding the crossed-box diagram (diagram B of Fig. |l[) to the kernel introduces no analytic structure into the 
problem beyond that discussed above, and therefore the Wick rotation may proceed exactly as for the ladder BSE. 
Indeed, the only change which needs to be made to our discussion of the ladder BSE is that Vi now includes a piece 
from the crossed-box diagram. 

In the case of the (fP'a field theory the value of the crossed-box diagram may be calculated using the usual Feynman 
parameterizations [ p9| , and the result is an amplitude 

1 l-Q l-a-/3 

167r2A2 



m.<i.^s)^lda I dp I ^7 ^,(^^^^^|^,^^^^^^) , (3.27) 



00 
where: 



D{a, f3, j\q'^, q^;s) ^ ajsn + P{1 - a - f3 - 7)322 
+ a(l - a - P - 7)s33 -I- /37S44 
-I- al3u + 7(1 — a — f3 — j)t 

- {a + P)m^ - {I- a~ (3)f/, (3.28) 



with: 
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^ 2 



■S33 = + 



10? 



90 j 



' ^2 



t = {qo- Qof - - + 2qq'q ■ q\ 
u^iqo + q'of - 9' - - 2qq'q ■ q' 



(3.29) 

(3.30) 

(3.31) 
(3.32) 



The amplitude X can thus be calculated numerically for any value of the parameters ^/s, q^, q' , qo, q, q ■ q' ■ There 
are no complications for ^/s < 2m + /z, since in that region the function D is negative for all relative energies for which 
the kernel must be evaluated. On t he ot her hand, above y/s = 2m + /x, singularities in the integrand are possible. 

The result of the integrations in (3.27) can be partial-wave expanded according to 



Xeiqo,q',qo,q;s) = J dxX{q'^,q^;s)PpX3 



(3.33) 



where x — q' ■ q, using another numerical integration. To obtain a phase shift of the same accuracy as in the ladder 
case 4 quadratures in each of the Feynman parameter integrations and 3 quadratures in the partial-wave projection 
are required. There is no need to change any of the quadrature numbers from the ladder solution. 

In the case fi = m our results agree with those of Levine and Wright |3^. In the case of interest to us, /i = 0.15, 
A — 0.13, the crossed-box diagram has a notable contribution. This is illustrated in Fig. ^ where we plot the phase 
shifts when including both dressing of the propagator and the crossed-box diagram (dotted line). For comparison 
the phase shifts resulting from the solution of the BSE in the ladder approximation with dressed (j) propagators are 
also included (dashed line). Comparing the results in Figs. ^ and ^ demonstrates very clearly the relative importance 
of the crossed-box diagram and the dressing of the (j) propagators. The contribution of the crossed-box diagram is 
significant in this case despite the small coupling used in the calculation. This suggests that the solution of the 
4>(t> — acfxj) equations will be substantially different from the ladder result — even below the cr-production threshold. 



IV. NUMERICAL SOLUTION OF THE <^<^ - a(j}<j) EQUATIONS 



The main motivation for examining the numerical solution of the (j)(f) — a(j)(j) equations is to develop the numerical 
approximations needed for the solution of the NN — ttNN equations, where spin and isospin are an additional 
complication. However, solving the (fxj) — acjxp equations will also allow a comparison of their solution with the results 
of the BSE with dressed (j) propagators, and either one-cr exchange or one-u exchange plus crossed-box diagram as its 
kernel. The analysis of Sec. || shows that the solution of the BSE is not going to be identical to that of the coupled- 
channels problem. However, a comparison may give us some insight into the relative content of the two equations and 
a measure of how important this difference is. This might tip the balance in favour of one or the other approach when 
considering NN scattering above the threshold for pion production. Even below the pion-production threshold, recent 
nucleon-nucleon interactions based on meson exchange have included the crossed pion exchange within the framework 
of two-body equations rather than in a coupled-channels approach |5iy52| ] . The results of this section might help us 
determine which of these two methods is best suited for including such crossed diagrams. 

In order to simplify matters as far as possible we have assumed that (f>a scattering is dominated by the pole diagram, 
i.e. t'^^j^ = 0. This yielded the two coupled integral equations (2.1C) and (2.11) with the potential V of Eq. (2.14). If 
we also assume that the (f) — 4> interaction is given by one-cr exchange and this in turn is represented by an S'-wave 
rank-one separable potential, then the input 0—0 amplitudei^ in Eq. ( 2.1C ) can be written as in Eq. ( ^.20| ). We 
may therefore recast the coupled equations, Eqs. ( 2.10| ) and ( 2.11 ), as 



7d0 



1 



Zdd> did2Td, 



(4.1) 
(4.2) 



where (j) refers to the (j)^) channel, while d refers to the a{(j)(p)d channel. Here Td is the propagator for a correlated 
(fxj) system. In writing the above coupled-channel equations we have taken into cons ideration the fact that the is 
a boson, and the equations need to be symmetrized. Although Eqs. (4.1) and (4.2) can be combined into a single 
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BSE for T^^^ we fiave chosen to solve the equations as a coupled-channels problem — in anticipation of the fact that 

if t^^^ 7^ the coupled equations cannot be reduced to a single BSE. 

In what follows all vectors are four-vectors unless otherwise stated. We define the relative four-momentum, q, of 
the d — a state of total four-momentum, P, via 



Pd = VdP + q 



where: 



Vd = 



Pa = VaP - q 



2s 



2s 



(4.3) 



(4.4) 



This choice for the relative momentum is made to ensure that when both the d and a particles are on-mass-shell 
go = 0. This in turn guarantees that the pinching of the q'^ integration contour in the BSE for {(f)(j))d(T scattering will 
occur at the point ~ 0, q " = q, where q is the magnitude of the on-shell three-momentum. The four- vector q which 
is defined by (4.3) and (4.4) is known as the Wightman-Garding relative momentum p3| . 
With this choice of kinematics we can write the input amplitudes Z as 



Zd4,{Pl,P2,Pl,P2) 
Z4,d{p'l,P2,Pl,P2) 



-z(2^)V4)(p;-fp^ 

Zd4,{Pl^P2,P'l,P'2) ■ 



- Pi - P2) [VAq - q') + Va{q + q') - X[q' , q; P) - X{q', -g; P)] 

- PI - P2)2 [Vd^q', q; P) + VdM, ~q; P)] , 



In Eqs. (|4.5|) and (4.6), the potentials Va for cr-exchange, and Vdtf, for 0-exchange, are given by 



Va{k) = 
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The 



Vd^{q',q;P)=g 



iSd{{i^dP + q?) 
{vdP + qY - ml ' 

We now define 

f^^ = (271)4^(4) + _ _ p2) (1 + P12) 



1 



fc2 — /i^ ' 
and a propagators are 

Td{q\P) - 



q? 



^p- 
2 



q 
2 



da{q\P) = 



yaP-q? 



/i- 



(4.5) 
(4.6) 
(4.7) 



(4.8) 



(4.9) 



fd^^-i{2nf5'^^\p[+p'2-pi-P2)Td^{l + Pi2) , (4.10) 



where P12 is the permutation operator for the two 0's and the factor of 1 + P12 ensures that scattering takes place in 
symmetric states only. In particular, no scattering can take place in a state with odd angular momentum. By contrast, 
if acting on a state with positive parity 1 -I- P12 merely gives a factor of two. The equations for these unsymmetrized 
amplitudes T^^ and Tdcf, are then 



T4,4,{q',q;P) 



+ i 



V^4q',q-P) + i 
d^q" 



'4„ll 



d'^q' 



(27r)4 



Td4q',q;P)^Vf^{q',q;P)+i 

where G4, is defi ned in Eq. ( |3.7| ), and Gd = 
Eqs. (4.11) and (4.12) are given by 



V^^iq', q"; P) G^{q"-P) T^^{q" , q; P) 

(4.11) 
(4.12) 

d system. The 'potentials' in 



(27r)4 

V^dW. q"\ P) Gd{q"; P) Td^iq", q; P) 
d^q" 



(2^)4 ^^^^"^'^ ^) ^^(9"; ^') T^^{q". q; P) , 



-Tddcr is the propagator for the a 



V^^{q'. q; P) = Vaiq' ^ q) ~ X{q' , q- P) , 

V^diq'. q; p) = V^diq', g; P) + V^d{~q', g; p) , Kf^(<?', g; P) = Vd^q', q; P) + Vd^{q\ -q; P) . 



(4.13) 
(4.14) 



If we remove the coupled-channel term Vd^ and the double -coun ting subtraction X the first equation here reduces to 
the Bethe-Salpeter equation solved for 00 scattering in Sec. [II A, thus justifying our use of the distinguishable-particle 
BSE to determine the 00 00 amplitude in the allowed two-body channels. 

Having derived these coupled integral equations they must be recast into a form convenient for computation. This 
requires three steps. Firstly, a partial-wave expansion must be performed in order to reduce the dimensionality of the 
equations from four to two. Secondly, any pinches causing difficulties in the energy range of interest must be removed. 
(Note that such pinches are not removed by Wick rotation.) Once this is done Wick rotation can be performed, thus 
moving the region of integration away from the poles of the propagators and exchange diagrams. 
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A. Partial-wave expansion 



The partial- wave expansion is performed exactly as defined in Sec . [II A. The amplitudes T^^ and T^^ and the 
"pot ential s" V^, X, Vd^ and V^d are expanded according to Eq. ( |3.5|) . This produces the results given in Eq. (3.g) 
and (3.33) for VJj and X. For Vd^, partial fractions may be used to split the integration, leading to: 



%/A 1 

D = 2q'q{jd^ - 7/3) 



ld4, 



2q'q 



113 



90 



2 / 



With this result in hand it is easy to compute V^d-i, since 

Vci,d-A%^ 9' 90, q; s) = Vd4,-eiqo, q, q'o,q'; s) 



(4.15) 
(4.16) 

(4.17) 



(4.18) 



The use of these partial- wave expansions leads to a set of coupled two-dimensional integral equations for the half-off- 
shell partial-wave amplitudes T^,p-i and Td^-/. These equations must then have all pinches removed. 



B. Kowalski-Noyes reduction 

Although the main interest in a coupled-channels approach is to investigate scattering in the energy region where 
three-body unitarity plays a central role, at this stage we are going to restrict our analysis to the energy region below 
the three-body unitarity cut. This is partly because in this energy domain we can use the Kowalski-Noyes approach, 
discussed in the last section, to remove the two-body unitarity cut. 

Since the coupled-channel equations do not couple amplitudes of different angular momentum, we now simplify the 
notation by dropping the angular momentum label on all amplitudes. 

The Kowalski-Noyes reduction is carried out for the two amplitudes in question by making the factorizations 

T^M,q',0,q;s) = UM^q')T^;{s) , Td4qiq\0,q;s) = fdM^q')T;;{s) , (4.19) 

where the dependence of the Kowalski-Noyes half-off-shell functions on s has been suppressed. 

Using these factorizations the partial-wave-expanded coupled-channel equations for the KN half-off-shell functions 
become 



oo 



yon 

— LXJ U 

OO oo 

-i J dq'^ I dq" %d{q'„ q\ q'^, q"; s) Gd{q'^, q"; s) fdM. q") , (4-20) 

-oo 
S t I I - \ °° °° 



fd^qi q') = ^ yon - » / dq'^ / dq" VdM, q', ql gM. 5) fM. q" 



-00 

00 00 



where: 



-I J dq'^ J dq" Vddiq'o,q', 9o , Gdiq'^, q"; s) fdM^ q") . (4-21) 

-oo 



t7 / / , \ -i, t ' > ^ y<t><t>iq'o^q'^0,q;s)V^^{0,q,qo,q;s) 

V4,4,{qo, q , go, 9; s) = V^^q^, q , 90, 9; s) , (4.22) 

,y f , \ ^,S( , > N ^00(90, g', 0, g; s) T//^(0, g, go, g; s) 

V4,d[qo,q, go, g; s) = F^d(go, g , go, g; s) tt^ , (4.23) 
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Vd0Wo) g. go, g; s) = V"rf0(go> g , go, s) ^75;;^ , (4.24) 

.yL, , ^ ^/0(go,g',O,g;s)"K^^d(O,9,go,g;s) 

V^idlgo, g, go, g; s) = . (4.25) 

It is clear from the structure of these Vs that the pinch in G^{qQ, q"; s) is now always ameliorated by a zero in V. 
In order to obtain phase shifts it is necessary to calculate the on-shell matrix, T^^[s), via 

= .^-tT^I-tI ^ ' (4.26) 
1 + 1/0(5) + i/d(s) 

with: 

00 00 

I^{s) = j dq'^ J dq" V^^{0, q, q'^, q"; s) G^iq'^, q"; s) UM^ g") , (4-27) 

-oo 

oo oo 

I4s) = J dq'^ I dq" V^diO, q, q'^, q"; s) G^q'^, q"; s) f^M^ g") ■ (4-28) 

-oo 



C. Wick rotation 

By writing coupled integral equations for the KN functions we have overcome the problem of the pinching between 
the poles of the (j) propagator. However, to convert the equations to matrix equations we have to avoid the other 
singularities of the kernel. This can be achieved by performing a Wick rotation in the integrals. Such a rotation 
requires that we take into consideration the singularities of both the kernel and the amplitude in the first and third 
quadrant of the (jg-plane. These singularities have been examined in detail p8| with the result that for 2m < -^/s < 
+ /i the only analytic structure which may cause problems during the rotation are the two poles from the Green's 
bmction G^. 

The final equations after Wick rotation can thus be written as 

oo oo 

■ / i\ ^0(*go, g , g, I I J II I J II -CrS I ■ I I ■ II ii\ I ■ II II \ £ I ■ II ii\ 

— + / c'go / dq V^^{iqo,q ,tqo,q )G^{iqo,q ; s) J ^^{iq^ , q ) 

■^■^ 

II -TrS I- -I I -( ll\ II \ 



"7 ^'^ ^00(^go.g.-(g )'g ) 2^sw{q")E,{q") '^'^'^ ) 

OO 

dgo / dq" %d{iq'o. g', «go , g"; s) Gdiiq'^, q"; s) U^iiq'^, q") , (4.29) 



S I ■ I I - \ °° °° 
f r I i\ K!(/)(*go, g , 0, <Z, s) f -I II f J ii-TrSc I I ■ II II \ri {■ II II \i I- II ii\ 
/d0(«go,g) = T7^ ^ I ^go / dq Vd^[iqo,q ,iqo,q ;s)G^[iqQ,q ; s) f^^(iqo , q ) 



. J ^g-^.-.(^go,g-,^(g"),g") ,^^(^.)^^(^.) g..(g") 

OO 

dq'^ f dq" VddiWo, g', *go , g"; s) Gd{tq'^, q"; s) fd4>{^q'^ , q") , (4.30) 



where in order to simplify the equations we have assumed undressed particles in taking the residue of G^, but the 
appropriate modification for dressed particles may easily be made. Also: 
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Vd0('?o, 90 , 9"; s) - g', g^,', g"; s) + VdM^i'^ , = 21/^^(9^,, g', q'^, q"- s) 

Note that g,pct,{q') = f<i)<p{w{q'),q') obeys the equation 



oc 00 



/ r\ Viptplwlq ),g,0,g;s) I 1 n I J II -TrS f - f i\ I ■ II II \ ri I ■ II ii \ t i ■ ii ii\ 
gMl)^ ^ ^ / dqo I dq V^^{w{q ),q ,tqQ,q ]s)G^[iq„,q ; s) f ^^[iq^ , q ) 



(4.31) 
(4.32) 



TT / dq"V^Aw{q'),q\w{q"),q"-s) 



2^sw{q")E^{q 



7K94"t>{<l ) 



00 00 



dq'i J dq" V^diwiq'), q\ iq^ q"; s) Gd{iq'^, q";s) U^{iq'^, q") . (4.33) 

-oo 

Wick rotation is also necessary to remove the poles in the integrals in Eq. ( 4.26| ) . Above we saw that 

14s) = 2a(s) + bis) , (4.34) 

where a{s) and b{s) are given by E qs. (3.15 ) and (|3.16| ), but with V^^ and /^^ replacing Fand /. As for Id, by Wick 
rotating the integration in Eq. (4.28) the result 



00 00 



ildis) = - J dq'i J dq" V^diO, q, iq'o, Gd{iq'^, q"; s)fd4iq'^, q"; s) 
-00 



(4.35) 



is found. In contrast to the BSE equation, where we had two coupled integral equations for the KN functions, here 
we have three coupled integral equations, with the a{(l)(f>)d channel being the source of the extra equation. 



D. Numerical results 



The equations (|^, ( p3o| ) and ( |43^ ) are in a form suitable for computation, having had almost all poles and cuts 
removed from the kernel. The only remaining pole is that at = 0, q" = q, due to the zero in V at that point being 
only first order, while the pole in G at that point is second order. However, as in the ladder BSE case, this pole may 
be dealt with by makin g a cha nge of integ ration variable and choosing appropriate quadratures. 

The equations ( 4.29 ), (OC) and ( 4.33 ) are sol ved f or the K o walsk i- Noyes half-o ff-sh ell functions and the results 
used to calculate phase shifts according to Eqs. ( 3.1C| ), ( 4.26 ), ( 3.15| ), ( |3.16 ) and ( 4.35 ). The equations are solved 
by discretizing the integrals using Gauss-Legendre quadratures and then applying matrix inversion. In order to 
guarantee phase shift accuracy to three significant figures 14 and 22 quadratures are necessary for the and q" 
(/></) intermediate-state integration, while 44 and 14 quadratures are necessary in the and g" integrations over the 



)d<^ intermediate-state. (Note that the integration is from -co to 00 in the ((/)</>)dcr case.) 8 quadratures are 
needed for the g" integration over [0, g]. 

As a first step in the application of the coupled-channel equations, we examine the contribution of the double- 
counting subtraction. In Fig. |l^ we present the phase shifts for </> — </> scattering with (solid line labeled C C — X) and 
without (dash-dot line labeled G G) this subtraction. Included for comparison are the phase shifts resulting from the 
solution of the BSE with (dotted line labeled Ladder-|-X) and without (dashed line labeled La dder) the crossed-box 
diagram. All four curves are calculated with the propagators dressed as described in Section [II C. 



Firstly, observe that the double-counting subtraction has a notable but not substantial effect on the phase shifts. 
This is a measure of the contribution of the crossed-box diagram, and the results are consistent with the difference 
between the ladder and ladder plus crossed-box diagram BSE calculations. Secondly, it is easily seen from Fig. |l^ 
that there is a considerable gap between the phase shifts given by the ladder plus crossed-box calculation (Ladder -|- 
X) and those produced by the coupled-channels calculation with the double-counting subtraction [C G — X). The 
difference between these two calculations is that G G — X contains the complete coupling to the (00)^(7 channel in its 
kernel, whereas Ladder -I- X includes only two Feynman graphs in its kernel. Thus we would expect some gap between 
the phase shifts predicted by these two calculations. However, at first sight it is odd that such a large difference could 
occur below the cr-production threshold, and with the small coupling strength, A = 0.13. 
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E. Effects of vertex dressing 



The explanation of the large difference could be due to the two diagrams (a) and (6) of Fig. g. These graphs are 
effectively included in the coupled-channels calculation, thus partially dressing the (fxjxr vertices. By contrast, the 
vertices used in the ladder and ladder plus crossed-box calculations were undressed. 

To gauge the size of this effect, we introduce vertex dressing in the ladder BSE and adjust the coupling strength A 
to get a binding for the </> — state that is comparable with the deuteron binding energy. We then use a separable 
approximation to this new ladder BSE amplitude as input to the coupled-channel equations. In this way diagrams (a) 
and (&) of Fig. ^ are included in both the BSE and coupled-channel equations. The effect of this on the coupled-channel 
equations is that some vertices in Fig. |^ are now bare vertices, and some are dressed vertices. 

We begin by seeking a self-consistent equation for the one-cr-exchange piece of the C C — X calculation. In other 
words, we must determine what form the acfxj) vertex function must have if it is to be the vertex in both the input 
(j)4> ladder interaction and the one-cr-cxchange piece of the (fxj) interaction which results from the coupled- channels 
approach. It is easily demonstrated that for this to be the case the vertex function / must obey 

f{pl,Pu,k,)=gb + igl / , , ■ -2 7- j— — —^f{p'^^k'^,p,^k'^,k,) 



(27r)4 k'^k"" - /i2 (p' - k')^{p' - k'Y - m2 •' (p - k')^{p - k'Y - m2 

(4.36) 



(This equation is given in diagrammatic form in Fig. |Tl|.) In general, equation (4.36) is difficult to solve. On the 



other hand, a reasonable approximation is obtained by arguing that the / in the integrand may be approximated by 
5, the dressed coupling constant. Such an approximation includes the one- loop dressing of the vertex exactly, and 
higher-loop dressings approximately. Its use may be justified by appeal to the same argument which validated our use 
of one-loop propagator dressing above. Higher-loop dressings make contributions to (f)(j)aa and higher-state unitarity. 
The one-loop result is therefore all we need to enforce 4>(j)a unitarity. Moreover, due to the small value of the coupling 
being used we might expect that any error in such an approximation is small. 

The integration over four-momentum k[, may now be done via the Feynman technique, and the result is 



npl,p.,k.)^g,-9X, I da I (4.37) 



D{a, f3\p'^,p^, k^) = /3(1 - « - /?)p>"^ + apk^k^ 

+a(l - a - /3)p^p'^ - (m^ - ^'^){a -I- /?) - M^ (4.38) 

where p^ is the initial nucleon four-momentum, p'^ the final nucleon four-momentum, /c^ the pion four-momentum, 

2 

and Ah = jf^, with the bare coupling constant. 

Defining / to be the value of the integral for p^ — p'^ = m? and k^ — fi^ (without the factors of g and gi, included) , 
suggests that the dressed and bare coupling constants are related by 

When the integral / is evaluated numerically we find / = —1.719. 

To include such dressing for the vertices in any of our previous BSE calculations only requires a change in the 
definition of the "potential" V. The additional vertex dressing does not interfere with the Wick rotation, since none 
of its cuts intrude into the first and third quadrants for ^/s < 2m + fi. The momentum-dependence of t he vertices 



may be calculated straightforwardly up to -^s = 2m + fj,, since at these energies the integrand in Eq. ( 4.37 ) is regular 
for all four- momenta used in the kernel of the integral equations for the KN functions. 

Numerical partial-wave expansion is used to evaluate the one-(T-exchangc interaction with dressed vertices. To 
obtain an accuracy of 3 significant figures in the phase shifts and 4 significant figures in the squares of the bound-state 
masses, 8 quadratures in this partial-wave integration and 4 quadratures in the Feynman integral is sufficient. The 
dressing of the 4> propagator is done as before but with a strength -y/AAb, rather than A, in order to simulate the effect 
of dressing one of the two vertices in the dressing loop. 

Firstly we consider the bound-state calculation. Tests show that the discretization of the integral equation can be 
carried out with the same number of quadratures without the accuracy diminishing. It is found that the bare coupling 
required for the scalar "deuteron" to have the desired mass is Xb — 0.0935, corresponding to A = 0.133. 
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To examine the importance of the momentum dependence of the vertex dressing, in Fig. ^ we compare the phase 
shifts with dressed propagators and vertices (sohd hne) with the a ladder calculation in which there is no vertex 
dressing and strength of the bare coupling is taken to be = A = 0.133 (dashed line). These results show clearly 
that the momentum dependence introduced as a result of the dressing only has a small effect on the final phase shifts. 
In fact, it reduces the phase shifts and to that extent increases the difference between these results and the results 
of the coupled-channels calculation reported in Fig. O. The small effect of the momentum dependence on the phase 



shifts provides justification for the approximation used in solving (4.37). 

We can now use the results of the BSE with dressed vertices and propagators as input to the coupled-channel 
equations. To do this, we need to fit the solution of the BSE in the ladder approximation with dressed vertices and 
propagators with a separable potential. The result of this fit is presented in Fig. ^ The best fit is achieved with 
j3 = 0.46760, which is not very different from the (3 found in the undressed-vertex case. 

This separable expansion of the ladder BSE with dressed vertices can now be used as the input to the coupled- 
channels calculation. However, if this is done additional double-counting is introduced, since the Feynman graph in 
Fig. [ij is double-counted, and so a further subtraction in V^^ is required. 

Other complications arise because the vertices which produce the coupling to the a{(j34>)d channel should be dressed. 
However, there is no (t>a non-pole t-matrix, and therefore no way for these vertices to be dressed in the coupled-channels 
theory derived above. To cover this deficiency of our calculation we arbitrarily replace gi, by g everywhere in V^d and 

The problem with this ad hoc solution is that the ((/)0)^cr coupling vertices which lead to vertex dressing for the 
one-a-exchange 4>(f> interaction are now also dressed. Therefore corrections are made to the interaction V^^, in order 
to remove all g's which lead to vertex dressing and replace them with (/f,'s once again. 

Note: dressed vertices in the crossed-box diagram, the (j) propagator dressing loop, and the "true" coupling to the 
{(j)(j))d'J channel all have no momentum-dependence — they merely contain the coupling constant g rather than gf,. 
We think this is a reasonable approximation to the full result, since the ladder results show that, for the values of 
the coupling of interest here, it is the overall strength of the coupling which is the major effect, not its momentum- 
dependence. 

Once these corrections to V^^, V^d, and Vdcj, are made the calculation proceeds exactly as above. The number 
of quadratures required for three significant figure accuracy remains unchanged. The results of these calculations 
are shown in Fig. |l^. Coupled-channels calculations both with and without the crossed-box diagram subtraction 
are compared to the ladder calculation with dressed vertices and the same calculation with the crossed-box diagram 
added. Now that graphs like (a) and (6) of Fig. |[ but with bare vertices, are included in the kernel of the BSE, the 
C C — X and Ladder -|- X calculations contain the same second and fourth-order Feynman graphs. Therefore, it is 
the two middle curves of Fig. ^ which should be compared. These two curves are much closer together than was the 
case in Fig. |l^, thus indicating that the discrepant vertex dressing was the main reason for the large gap between the 
same curves in that plot. 

So, we conclude that if the vertex dressing is done approximately consistently in the Ladder -I- X and C C — X 
calculations then the resultant phase shift curves lie very close together. The double-counting subtraction is crucial 
to our obtaining this agreement. 

V. CONCLUSION 

In this paper a coupled-channels formulation of (jxf) scattering in a (jP'a field theory was investigated. In particular, 
we found that the ipcj) — ucjxf) equations can be solved using an extension of the method developed by Levine et al p^ ] 
for the ladder Bethe-Salpeter equation (BSE) . This may be done in a moderate amount of computer time on a SUN 
workstation. Since these equations are just the NN — nNN equations of Ref. adapted to a scalar field theory with 
i^^j = this indicates that it is feasible to solve the four-dimensional covariant NN — t:NN equations numerically. It 
also yields phase shifts for (fxj) scattering up to the {(j)(j))dO' threshold, thereby showing that the (jxp — acjxj} equations can 
be used to effectively include in the BSE kernel the infinitely many diagrams involving one explicit a. So, it allows 
a comparison between such a coupled-channels (CC) approach to (jicj) scattering and descriptions using a single BSE, 
which allow only a finite number of Feynman graphs to be included in the kernel. 

As discussed at length in Refs. [^,^,^, in most previous four-dimensional NN — ttNN equations some diagrams 



are included more than once in the kernel of the coupled-channel equations. In the first calculation of Sec. IVD 
the only double-counting removal required is that of the crossed-box diagram. Its subtraction is found to make a 
significant difference to the phase shifts obtained from the coupled-channels calculation. 

However, the gap between the phase shifts produced by this, C C — X, calculation and those obtained when the 
crossed-box diagram is added to the kernel of the ladder BSE (Ladder + X) is surprisingly large. Closer examination 
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reveals that the 4>4>'^ vertices in the two approaches are not dressed in the same way. The coupUng to the (00)(jcr 
channel in the C C formulation introduces some dressing of the vertices. If a direct comparison is to be made this 
dressing must be included in the ordinary BSE calculations. Once this consistent dressing is implemented a striking 
decrease in the gap between the C C — X and Ladder + X calculations occurs. 

It is then seen that the "true" coupling to the {(j)(f>)df7 channel makes little difference to the phase shifts. In other 
words, if all of the second and fourth order diagrams which are effectively included in the C C — X calculation's kernel 
are summed in the kernel of the single BSE then the results of the two calculations are very similar. It is worth noting 
that the agreement between the single and coupled BSE approaches is destroyed by the omission of any fourth-order 
diagram from the Ladder + X calculation, and by the over-counting of the fourth-order crossed-box graph in the C C 
calculation. Whether this close agreement persists above the {(f>4>)dU threshold remains to be investigated. 

It is only because the "true" {4>4))d(J channel is relatively unimportant that we can claim to have summed the field 
theory correctly. If this were not the case our use of the ladder BSE as the basis for the input (jxp interaction would be 
open to criticism. Strictly speaking in order to "bootstrap" the theory up we should construct a parameterization of 
the (jjcj) amplitude resulting from our C C — X calculation and use this as input to a new C C — X calculation, repeating 
this process until convergence is obtained. However, if this procedure were implemented additional double-counting 
would be introduced into the equations. This would have to be explicitly removed. As the "true" ((j)(f>)d(7 coupling 
appears fairly small in our calculation, we do not pursue such a bootstrapping procedure here. 

In this paper we have only calculated cfxp scattering up to the cr- "deuteron" threshold. In modifying the methods 
of this paper for work above ^/s — nid + fJ- three issues arise. Firstly, once this energy is reached the a and d 
poles pinch each other, thus generating the elastic {(t)4>)d<J threshold. This pinch and threshold may be removed by 
modifying the Kowalski-Noyes method so that it applies to coupled-channels problems. Secondly, the validity and 
efficacy of Wick rotation must be examined. The discussion of analytic structure in Appendix G of Ref. |28| shows 
that, for rud + A* < i/i < 2m -|- 2/z, the only cut which threatens Wick rotation is the one which represents the process 
{4"t')d However, even if y/s is large enough for this "deuteron" break-up cut to intrude into the third quadrant 

the cut may still be avoided by Wick rotating about a point on the negative Qq axis in {(j)(l))d<J intermediate states. 
Provided this point is suitably chosen, the only structure which can obstruct such a rotation is the cr-propagator 
pole. Hence an auxiliary equation for the Kowalski-Noyes functions at q'^ = v^^fs — E^{q") must be written when 
the Wick rotation is performed. Once the rotation is completed the only new analytic structure in the kernel occurs 
in that part of the kernel corresponding to the transition {(j)(l))d + c (jj + (j). In the auxiliary equation for /^^ at 
iij{q') that piece of the kernel acquires logarithmic branch-points, due to the possibility of real a production via the 
diagram shown in Fig. |l6|. Thus, Wick rotation is permitted and succeeds in eliminating most of the troublesome 
analytic structure from the kernel. Finally, the increase in available energy leads to complications in the calculation 
of Feynman integrals. In both the crossed-box graph and the vertex loop some of the intermediate-state particles can 
propagate on-shell once y/s reaches 2m + fi. So, singularities occur in the integrands of the Feynman integrals and 
the integrals can no longer be done in a straightforward way. 

Another point for future investigation is the inclusion of a non-pole (pa interaction in the theon^ This might be 
done by solving the ladder BSE for unequal mass particles using the driving term shown in Fig. |l^. The resultant 
amplitude could then be parameterized by a separable interaction and used as input to the (jycj) — acfxj) equations. 
This obviously increases the amount of physics summed in the theory, but also has three other effects. Firstly, by 
introducing an additional channel, (0cr) -t- 0, the computer time required is raised. Secondly, the vertex dressing 
problem no longer needs to be remedied "by hand" since all vertex dressing will be done in a consistent way. Finally, 
further double-counting corrections need to be introduced, as discussed in Section II. However, since all of these 
further double-counting corrections involve the removal of the crossed diagram from the (pa interaction this third 
point should not complicate matters greatly. 

Few changes will be needed in order to apply the methods of this paper to the NN — ttNN system. The first step 
in such a calculation will be to perform the spin and isospin algebra involved in the partial-wave decomposition of the 
four-dimensional NN — ttNN equations. Once this is done a set of coupled partial-wave expanded equations of the 
form solved in this paper is obtained. Questions remain about the validity of Wick rotation in the NN — ttNN case, 
since some of the propagators now are proportional to i, rather than the of the scalar case. But, the presence 
of form factors in our equations will ensure that the integrands go to zero fast enough to validate the use of Wick 
rotation. Thus, once a set of coupled equations for the NN — ttNN system has been obtained the approach described 
here may be applied. The computing time required for such a calculation will obviously be much longer than for those 
in the scalar field theory. It should also be pointed out that if results for the NN — irNN system above the two-pion 
threshold are needed a method which is more sophisticated than our "naive" Wick rotation should be used. 

However, for calculations up to the second production threshold the work of this paper shows that the numerical 
solution of the coupled field-theoretic equations derived in Ref. is entirely possible. 
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FIG. 1. The subtractions in V required to avoid double counting. 



FIG. 2. The lowest-order Feynman diagrams that contribute to Veff if the input — <^ amplitude is the solution to the ladder 
BSE. 



FIG. 3. A pictorial representation of the ladder Bethe-Salpeter equation as given in Eq. (3.2) 



FIG. 4. Coupling constant versus bound-state position for the ladder BSE with m = 1 and /i = 0.15. 
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FIG. 5. Comparison of the separable approximation (dashed hne) to the ladder calculation with undressed particles (solid 
line) for jj, = 0.15 and A = 0.13. 



FIG. 6. The Schwinger-Dyson equation for the 0-particIe propagator. Heavy lines represent dressed propagators, and the 
lighter lines represent bare propagators. 



FIG. 7. Real part of the S'-wave phase shifts in the ladder approximation, with (solid curve) and without (dashed curve) 
one-loop dressing for the case A = 0.13. 



FIG. 8. Imaginary part of the /S'-wave phase shifts in the ladder approximation, with (solid curve) and without (dashed 
curve) one- loop dressing for the case A = 0.13. 



FIG. 9. Real part of the iS-wave phase shifts with (dotted curve) and without (dashed curve) the crossed-box diagram, for 
the case of dressed particles with A = 0.13. 
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FIG. 10. Comparison of four different ways of calculating the real part of the S-wave phase shifts. The dashed curve is the 
result of the ladder calculation, the dotted curve includes the crossed-box diagram, the dot-dashed curve is the straightforward 
coupled-channels result, and the solid curve is the coupled-channels result with the crossed-box diagram subtracted so as to 
remove double-counting. 



FIG. 11. The non- linear equation for the vertex function /. 



FIG. 12. The real part of the iS-wave phase shifts for the case A;, = 0.0935 with vertex dressing included (solid line), compared 
with the case A = 0.133 with no vertex dressing (dashed line). 



FIG. 13. The real part of the 5-wave phase shifts for the case A;, = 0.0935 with vertex dressing included (solid line), compared 
with a separable approximation to this result (dashed line). 



FIG. 14. A Feynman graph which is double-counted in the coupled-channels calculation if an input 4> — 4> amplitude with 
dressed vertices is used. 
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FIG. 15. Comparison of single-BSE approach, with vertices approximately dressed to all loop orders, against coupled-channels 
approach. Legend as in Fig. HO. 



FIG. 16. The diagram which, above a-production threshold, leads to logarithmic branch-points in the kernel of the 
Wick-rotated coupled integral equations. 



FIG. 17. A possible driving term for the ladder BSE for (j><^ scattering. 
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